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Abstract. We show that for an eUiptic divisibihty sequence on a 
twist of the Fermat cubic, + — to, with m cube-free, all the 
0^ . terms beyond the first have a primitive divisor. 



in 



X 



1. Statement of Main Theorem 



^ ' Let C denote a twist of the Fermat cubic, 

C: U^ + V^ = mW^ (1) 

with m a non-zero rational number. If K denotes any field of char- 
acteristic zero, the set C{K) of projective i^-rational points satisfying 
J> I ([1]) forms an elliptic curve. With respect to the usual chord and tan- 

gent addition the set C{K) forms a group . The identity of the group 
is (-1,1,0) and the inverse of the point {U,V,W) is {V,U,W). Let 
\ R G C(Q) denote a non-torsion rational point. Write 

nR={Un,Vn,Wn), f/„, K, G Z 

in lowest form with gcd(t/„, 14,, W^n) = 1- This paper is devoted to 
proving the following theorem. 



O 

^ \ Theorem 1.1. Let C denote the elliptic curve in ([T]) with m G Z 

\ assumed to he cube-free. Let W = (Wn) denote the sequence obtained 

as above from R G C{Q), a non-torsion rational point. For all n > 1, 
the term Wn has a primitive divisor. 



The term primitive divisor is defined in the following way. 

Definition 1.2. Let {An) denote a sequence with integer terms. We 
say an integer c? > 1 is a primitive divisor of An if 

(a) d I An and 

(b) gcd{d, Am) = 1 for all non-zero terms A„i with m < n. 
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The sequence W = (Wn) is a divisibility sequence, which means that, 
for all m, n G N, 

m I n implies Wm \ Wn- (2) 

In line with recent developments [ni[l2l[T3l[Tl[l5l[T6l[25l[26l[27] 

we define the sequence W = (Wn) to be an Elliptic Divisibility Se- 
quence. Admittedly this stretches the definition originally used by 
Morgan Ward [29j but we believe it is a reasonable name for a di- 
visibility sequence that arises from an elliptic curve. 

In the sequel, it will often be convenient to work with the affine curve 

+ = m (3) 

and we will also refer to this curve as C. Properties of points such as 
being integral will generally refer to the affine curve. Thus, an integral 
point on C is a pair of integers {u, v) satisfying ([3]). Theorem II. II is the 
best possible in the sense that if R is an integral point, then Wi has 
no non-trivial divisors and so (by our definition) no primitive divisors. 
Thus, the identity 

(l + t)3 + (l-t)=^ = (6^2 + 2), 

for any integer t > 1, gives rise to an elliptic divisibility sequence for 
which Theorem 1 1.1 1 is best possible, whenever 6t^+2 is cube-free. There 
are infinitely many t such that + 2 is cube-free by a result of Erdos 

i- 

If one removes the condition that m be cube-free, then it is easy to 
construct counter-examples to Theorem 11.11 by clearing denominators. 
However, any elliptic curve in the form ([1]) can be transformed into 
one with m a cube-free integer by a simple scaling. Thus, Theorem 11.11 
does give complete information, taking account of the transformation. 

Note that Theorem 11.11 has an immediate application to the study 
of integral points on elliptic curves in the form ([1]). In [23], Silverman 
showed that, when m is cube-free, there exists an absolute constant k 
such that ([T]) has at most ft;i+i'ank(c/Q) integral points. In particular, the 
number of integral solutions to is bounded uniformly if we restrict 
attention to curves C with rank not exceeding a given bound. Note that 
C(Q) is torsion- free when its rank is positive. Thus if rank(C/Q) = 1, 
then C(Q) consists solely of the multiples of a single point R, say. 
By Theorem 11.11 the denominator of nR, for each n ^ ±1, 0, has a 
primitive divisor. In particular, none of these denominators is 1. 

Corollary 1.3. If m ^ Z is cube-free, and rank(C/Q) = 1, then ([3]) 
has at most 2 integral solutions - either of which generates C(Q). 
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Corollary 11.31 is not overly surprising. The method of proof in [23] 
gives very strong bounds on the number of integral points on C as the 
diophantine approximation involved in this case is trivial. Corollary 1 1.31 
is noted because it is sharp and completely qualitative. For quantitative 
results of similar strength for other elliptic curves, see |15j . 

In the next section, we will set Theorem II. II in its proper context, as 
well as outline the structure of the proof. The proof occupies the rest 
of the paper. 

2. Primitive prime divisors 

Let M = (M„) denote the Mersenne sequence, whose ?7.th term is 
M„ = 2" — 1. No proof is known that M contains infinitely many 
prime terms. The concept of a primitive divisor was introduced as a 
way of showing that new primes are produced by the terms of M, but 
in a less restrictive sense. In 1886 Bang [2j showed that if a is any 
fixed integer with a > 1 then the sequence with nth term a" — 1 has 
a primitive divisor for any index n > 6. This is a sharp result because 
the term Mq = 63 = MI.M3 does not have a primitive divisor. Bang's 
theorem is remarkable because the number 6 is uniform across all a and 
it is small. Indeed, it is not hard to show that a = 2 is the only example 
realizing this bound. Bang's Theorem was incredibly influential. 

In 1892 Zsigmondy [201 obtained the generalization that for any rel- 
atively prime a and b with a > 6 > 0, the terms of the sequence 
An = — 6"" all have primitive divisors if n > 6. This lovely result was 
re-discovered several times in the early 20th century and it has turned 
out to be quite applicable. See [20] and the references therein where 
applications to Group Theory are discussed. For example, for fixed q 
a prime power, let denote the finite field with q elements. Zsig- 
mondy's Theorem applied to the explicit formula for the order of the 
group GLn{¥q) shows this order has a primitive divisor for all large n. 
Thus Sylow's Theorem can be invoked to deduce information about 
the structure of the group. 

Definition 2.1. Let A = (An)^^^ be an integer sequence. Define 

Z{A) = maxjn | An does not have a primitive divisor} 

if this set is finite, and Z{A) = cxd if not. The number Z{A) will be 
called the Zsigmondy bound for A. 

Thus Bang's Theorem may be stated: Z{M) = 6, while Zsigmondy's 
Theorem may be stated: Z{A) < 6, where = a" — 6" as above. 
Following Zsigmondy's Theorem, the next major theoretical advance 
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was made by Carmichael. Let u and v denote conjugate quadratic 
integers. Consider the integer Lucas sequence U defined by 

Un = {u^-vn/{u-v). 

For example, the Fibonacci sequence F = (Fn) is a Lucas sequence. 
Carmichael [5] showed that if u and v are real then Z{U) < 12. This 
too is a sharp result because F12 does not have a primitive divisor. 
The general case was settled by Bilu, Hanrot and Voutier [3]. They 
proved that Z{U) < 30 using a powerful cocktail of methods including 
start of the art bounds from Diophantine analysis as well as massive 
computations to deal with special cases. Again this is a sharp result 
as the sequence generated by the roots of the polynomial — x + 2 
illustrates. The paper [3j gives details about the long journey from 
Carmichael's result to the general case. 

2.1. Elliptic Curves. Now let E denote an elliptic curve in Weier- 
strass form, 

y'^ + aixy + a^y = x^ + 02^^ + a^x + ag 

with ai,...,aQ G Z. The shape of the defining equation forces the 
denominator of x{Q) to be an integer square, for any Q G E{Q). Let 
Q G E{Q) denote a non-torsion point. For every n G N write 

x{nQ) = An/Bl 

in lowest terms. The sequence B = {Bn) is an elliptic divisibility 
sequence associated to Q and E. Silverman |2H] obtained a primitive 
divisor theorem for elliptic divisibility sequences arising from curves 
in Weierstrass form. It seems likely that a uniform version of this 
theorem holds for curves in global minimal form. In other words, if 
B = (Bn) arises from a rational point on a curve in global minimal 
form then Z{B) < Nq, where A'^o is independent of E and Q. The 
proof of Silverman's Theorem suggests that the Zsigmondy bound is 
higher for sequences generated by rational points with small global 
canonical height (therefore it is significant that all such heights are 
uniformly bounded away from zero). The following example appeared 
in [H]. 

Example 2.2. Let B denote the elliptic divisibility sequence generated 
by the point Q = (7107, —602054) on the elliptic curve 

y^ + xy + y = x^ + x'^ ~ 125615x + 61201397 

Computations suggest that Z{B) = 39 and no higher value of Z{B) is 
known for an elliptic divisibility sequence coming from a Weierstrass 
curve in minimal form. 
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The curve in Example 12.21 was taken from a list of small height 
points maintained by Noam Elkies The curves in Elkies' table 
are not generally in minimal form but the curve in Example 12.21 has 
been rendered in minimal form in order to estimate Z{B). 

Computations with congruent number curves suggest the Zsigmondy 
bound is generally very small. A uniform Zsigmondy bound appears 
in [12] and [H] for an infinite class of sequences arising from congruent 
number curves. Specifically, let T > 5 denote a square-free integer and 
let Q denote a non-torsion rational point on the curve 

2 3 T-i2 

y = X — I X. 

In [12] it was shown that, if x{Q) < or x{Q) is a rational square, then 
Z{B) < 21. In [14j this result was improved by reducing the Zsigmondy 
bound to 2 (and allowing any of x{Q) or x{Q) ± T to be a square), 
a bound witnessed by an infinite family of sequences. Provided the 
rank of the curve is positive, there will always be points satisfying the 
hypotheses stated: this is because, if x{Q) > then x{Q + [0,0]) < 
and, for any rational point Q, 

x{2Q) and x{2Q)±T 

are all rational squares. 

2.2. Comparisons with the classical theory. There are notable 
similarities between the results for elliptic divisibility sequences and the 
classical results described earlier. Both give a uniform bound across 
infinitely many sequences which is best possible, both rely upon lower 
bound on heights and both reduce the problem to solving a finite num- 
ber of Thue-Mahler equations. Also, the uniformity result for Lucas 
sequences relies on good bounds from transcendence theory together 
with the fact that the answer to Lehmer's problem is known for qua- 
dratic integers. A uniform result for elliptic divisibility sequences in 
general would appear to require better elliptic transcendence results 
than are currently known, together with a proof of Lang's Conjecture 
on lower bounds for heights of points on elliptic curves. 

In one respect, however, the arithmetic of these two classes of se- 
quences differs markedly. Bang's Theorem may well have arisen as 
part of an attempt to prove the Mersenne Prime Conjecture (which 
remains open). On the other hand the analogue of that conjecture is 
false for elliptic divisibility sequences on curves of the form ([1]), and 
heuristics, as well as proofs in special cases, indicate that it fails for 
those on Weierstrass curves too, see [10] and [TT] . 

In the section that follows we will reduce the given problem to one 
on the Weierstrass model of a curve birationally equivalent to the 
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curve ([T]). The method proceeds in a pincer movement, somewhat 
similar to that in the two papers [12] and [H]. These papers used a 
good lower bound for the canonical height of a rational point which 
were obtained in [1]. Here, our workhorse is the paper [l7j, although 
the height bounds are not stated or used in the same way as in [T7] . 
However, beside the similarities, there are many intriguing differences. 
Most notably, in this paper we make heavy use of a numerator se- 
quence on a Mordell curve, see (jlj), for which we can prove a uniform 
Zsigmondy bound - see Theorem I2.3[ Remarkably, this sequence is not 
a divisibility sequence, one of the few known cases where a primitive 
divisor theorem can be proved for a sequence which lacks the divisibil- 
ity property. Also, remarkably, we have been unable to prove a uniform 
Zsigmondy bound for the corresponding denominator sequence on the 
Mordell curve. 

2.3. The structure of the proof. The proof of Theorem 11.11 relies 
upon upon two different techniques. The first one bounds Z(W) above 
by showing that the non-existence of a primitive divisor of Wn implies 
a certain divisibility relation (see fl24l) ) involving the terms An, which 
are defined in This relation leads to an inequality which is violated 
for all sufficiently large n. This part results in an upper bound for n 
of n < 14. 

The second step shows directly that for each of the indices n > 2 
not covered by the first part, Wn does have a primitive divisor, by 
reducing the checking to a finite number of Thue-Mahler equations. 
In the reduction, the elliptic division polynomials, which are elliptic 
analogues of the cyclotomic polynomials, play a starring role. 

The combination of techniques described here runs exactly parallel 
to those used in earlier primitive divisor theorems such as those of |21 
[3l[5l[30]. Note that it is essential to reduce the bound for Z(W) in the 
first step as low as possible, in order to keep to a minimum the number 
of Thue-Mahler equations which need to be solved in the second step: 
without adequate care, computationally infeasible problems result. 

The proof uses the well-known bi-rational equivalence of ([3]) with the 
Mordell curve 



E -.Y^ = X^- 432m2. 



(4) 



The map is given by 



36m -|- Y 



and V 



36m — Y 
OX 



u = 



6X 
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If -R G C(Q) corresponds to Q E E{Q) under the transformation, and 
we write 



52' 53 

n n 



then 



Wn QAnBn '''' Wn 6A„5„ ' ^ ' 

The proof of Theorem 11.11 exploits both the denominator sequence 
B = (Bn) in ([5]), as in [l2l[14], and the numerator sequence A = (An). 
The latter may have independent interest. 

Theorem 2.3. Let E denote the Mordell curve and suppose Q is 
a non-torsion point in E{Q). Let A = [An) denote the sequence as 
defined in Then Z{A) < 12. 

In principle, techniques analogous to those used in [H] and Section [5] 
can be used here to reduce the bound stated in Theorem 12.31 In prac- 
tice, however, the computations involved for some cases are beyond 
our current capabilities. We are uncertain about the supremum of the 
values Z{A) as m varies. Perhaps it occurs when m = 7: in this 
case Z{A) = 2 because A2 is a proper divisor of Ai, whereas all the 
terms ^3, A4, . . . have a primitive divisor. 

3. Local arithmetic 

Since we are interested in the prime divisors of the numerators and 
denominators of x{nQ), for Q G E{Q), it makes sense to address the 
local arithmetic of E. 

Standing assumption: throughout Section [3|, p will be a 
prime other than 2 or 3. 

We will analyze E{Qp), as in [2U Chapter VII], through the exact 
sequence 

^ Ei(Qp) ^ Eo(Qp) ^ Ens(Fp) ^ 0, (7) 

where i?o(Qp) is the subgroup of E{Qp) consisting of points with non- 
singular reduction modulo p, Ei{Qp) is the kernel of reduction modulo 
p, and i^ns(IFp) is the group of nonsingular points on the curve reduced 
modulo p. Throughout the paper, the terms nonsingular reduction and 
good reduction will be used synonymously, as will the terms singular 
reduction and bad reduction. 

For an arbitrary integer sequence A = (An), follow Morgan Ward's 
terminology [29j and define the rank of apparition of the prime p in 
the sequence A to be the least index n such that p \ An (the rank is 00 
if no such index exists). For the sequence B defined above, note that 
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p I -B„ is equivalent to OTdp{x{nQ)) < 0, in other words, nQ G Ei{Qp). 
The rank of apparition of p in the sequence B is, then, the order of 
Q in E{Qp)/ Ei{Qp), and p \ precisely when k is divisible by this 
order. In fact, the power to which p divides -B^n is entirely predictable 
once we know the power to which p divides B^- The following lemma 
is obtained in [21] by appeal to formal groups (and by an examination 
of division polynomials in [H]). 

Lemma 3.1. Let Q e ^i(Qp). Then 

oidp{x{kQ)) = oTdp{x{Q)) — 2oTdp{k). 

Remember that p > 3 is a standing assumption in this section. 
Lemma 13.11 is not generally true when p = 2 (although it is true when 
p = 3). 

Note that the commentary above about the rank of apparition leads 
to an entirely algebraic interpretation of primitive divisors in the se- 
quence B. To say that the term Bn fails to have a primitive divisor 
is to say that there is no prime p such that Q has order exactly n 
in E{Qp)/Ei{Qp). This quotient, for primes p of good reduction, is 
simply E{¥p). The rank of apparition of p in the sequence A may sim- 
ilarly be interpreted in terms of the local arithmetic of E, although the 
interpretation depends on whether E has good or bad reduction at p. 

If p is a prime of bad reduction for E, then p \ m (recalling again 
that p > 3 here). The curve E has a singularity at the point [0,0], 
and so p \ An precisely when nQ is singular modulo p, in other words, 
nQ has non-trivial image in the quotient E{Qp)/ Eq{Qp). Since the 
discriminant of E is A(i?) = — 2^^3^m'^, we see that OTdp{A{E)) is 
divisible by 4 for all p | m other than 3, and one can easily check 
(see Table 15.1] or use the addition formula) that 

E{Qp)/EoiQp) = Z/3Z. 

Thus if p appears in the sequence A at all, then it appears in precisely 
the terms A^ for which 3 \ k. The primes so occurring are, of course, 
exactly the primes of bad reduction dividing Ai. 

If, on the other hand, p \ m, then E{Qp) = -E'o(Qp) aiid E = E^s- 
Let Hp C E(¥p) be the subgroup generated by the two points in E(¥p) 
with x-coordinate 0. It is easy to verify that Hp = Z/3Z, and clearly 
p I An if and only if 

nQeHp\ {O} 

where Q is the image of Q in E{¥p). Thus the rank of apparition of 
p in A is the order of Q relative to Hp \ {O}, or the least n such that 
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nQ e Hp\ {O}. Note that as -A32m'^ = -3{12mf, we have 
Hp C E{¥p) if and only if p = 1 mod 3, 

that is, precisely if E has ordinary reduction at p. Super-singular 
primes cannot appear at all in the sequence A, in marked contrast 
with the situation for the sequence B, in which every prime eventually 
occurs. Of course, even p being a prime of ordinary reduction for E 
does not ensure that p has finite rank of apparition in A. 

The usefulness of Lemma 13.11 from our perspective, is that it allows 
one to obtain a strong bound on the size of a term i?„ failing to have 
a primitive divisor, as in fl2[ [HI [16]. The main goal of this section is 
to prove a similar result for the sequence A. 

Lemma 3.2. Suppose OTdp{x{Q)) > 0. Then, for any k coprime to 3, 
OYdp{x{kQ)) = oTdp{x{Q)) + ordp(/c). 

Proof. Much of what we need to prove our result has already been 
established. If p is a prime of good reduction for E, then the condition 
OTdp{x{Q)) > ensures that Q E Hp \ {O}. In particular, 3Q = O, 
and so 3Q e i?i(Qp). As 3 f fc, we have kQ E Hp\ {O} as well, and 
consequently 3kQ G Ei{Qp). Indeed, triplication on E follows the law 

x9(Q) + 293V(g)m^ + 2i^3V(g)m^ - 2'^3''m'^ 
9x2(g)(x3(g) -2633m2)2 ' 

and so we deduce (for primes p \ m) that 

ordp(x(g)) > ^ ordp(x(3g)) = -2 ordp(x(g)). 
Applying Lemma [3. H we obtain 

2 ordp(x(A;g)) = — oidp{x{3kQ)) 

= - ordp(x(3g)) + 2 ordp(fc) 
= 2ordp(x(g)) +2ordp(A;). 

The proof is somewhat more involved for the case where p is a prime 
of bad reduction. In this case we know that E{Qp)/ Eq{Qp) = Z/3Z, 
and ordp(a;(g)) > ensures that Q is non-trivial in this quotient. If 
we let 

EniQp) = {g G E{Qp) : ordp(x(g)) < -2n}, 
then for each n > 0, 

EniQp) /En+i{Qp) = ^/P^- (9) 

This agrees with our definitions above of Eq and Ei, and the fact that 
E has additive reduction modulo p. It is worth mentioning that (Q 



10 GRAHAM EVEREST, PATRICK INGRAM AND SHAUN STEVENS 

is essentially equivalent to Lemma 13.11 (although ([9]) only holds in the 
general case for n > 1). 

Returning to (jH]) we see that ordp(x(Q)) > implies 

2ordp(x(Q)) = -ordp(x(3g)) +2ordp(m). (10) 

If 3Q G EniQp) \ En+iiQp) and r = ordp(A;), then tells us that 

3kQ e En+r{Qp) \ K+r+l(Qp) SO, by ([ID]), 

2oTdp{x{kQ)) = — ordp(x(3A;(5)) + 2 ordp(m) 

= — ordp(x(3Q)) + 2 ordp(A;) + 2 ordp(m) 
= 2oTdp{x{Q)) + 2oTdp{k). 

□ 

4. The proofs that Z{A) < 12 and Z{W) < 14 

After gathering some preliminary results, we proceed with the proof 
of Theorem 12.31 The proof that Z{W) < 14, which advances along 
similar lines, follows. The step from Z{W) < 14 to Z{W) < 1 is taken 
in Sections [5] and [61 

4.1. Preliminaries. It will not generally be true that the defining 
equation ([4|) for E is in global minimal form. The following comes 
from [Ul Lemma 1]. 

Lemma 4.1. //9 | m write M = m/9. Then the global minimal form 
E* for E is 



(I) 


E* : 




= x^ - 


27m2 
4 




m 


and 9 \ m 


(11) 


E* : 


y^ + y 


= x^ - 


27m2 + 1 
4 




\ m 


and 9 \ m 


(III) 


E* : 




= x'- 


3M2 
4 




m 


and 9 m 


(IV) 


E* : 


y^ + y 


= x'^ - 


3M2 + 1 
4 


^/2i 


\ m 


and 9 m. 



The following explicit transformations render the curve E in minimal 
form: 

X = u^x and Y = u^y + t 
where [u,t] = [2,0], [2,4], [6,0], [6, 108] (respectively). 

In the sequel, these four possibilities will be referred to as Cases 

(i)-(iv). 
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Lemma 4.2. Write Q for a non-torsion point on E{Q), corresponding 
to Q* G E*{Q), where E* denotes the minimal model. Write, for all 
n> 1, x{nQ*) = an/b'^n and h = h{Q). Also, write M = m/9 if 9 \ m 
and M = m otherwise. Then 



--logM 
3 ^ 



- log 3 < hn^ 



<^log3. 



The proof follows immediately from [TTl Proposition 2]. Note the 
misprint in [17j (which has + signs on the left hand side). We are 
going to use this in the following form: 



hn^ - — log 3 - - log 



1 + 



54M2 



ma;; 



1 2 2 

- log an<hn + - log M 



< 2 log ^« 

- - log 3. 

2 ^ 



Lemma 4.3. Let P denote any non-torsion point in 

Hp) > ^logm- ^log2- ^log3 > ^log 

unless m = ±2 mod 9 and m has a prime factor congruent to 1 mod 6, 
in which case 



and (11) 
(12) 

'(Q). Then 
m-.0562 (13) 



h{P) > — logm 



1 

27 



log2 - ^ log3 > ^ logm - .1173 (14) 



Remark 4.4. The difference between the bounds in (fT3|) and f|T4l) 
might seem so slight as to be hardly worth mentioning. However the 
sieving allowed by the second bound greatly reduces the amount of 
manual checking in the sequel. Having said that, it will become clear 
in the following proof that further sieving is possible. The lemma as 
stated represents a compromise between further savings on the checking 
of values of m as against a more complicated version of Lemma 14.31 

Proof. The proof of Lemma 14.31 uses the analysis in [T7] . It begins 
by estimating a lower bound for h{kP), for A; = 2 or 3, then uses 
h{kP) = k'^h{P) to obtain the bound sought. The global height is 
bounded by estimating the local canonical height Ap at each place. 
Write 

h{P) = Y,K{P)- 

p<oo 

We write Q = kP and estimate Xp{Q) for each finite prime p. Sup- 
pose that Q G Eo{Qp), the non-singular part of the p-adic curve; then 
from [17j (2.5)] we get 



KiQ) = 2 logmax{l, \x{Q)\p} 



1 

12 



log I A* 



IP! 
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where A* denotes the discriminant of the minimal equation, so 

KiQ) > -^log|A*|p. (15) 

For the archimedean valuation, assume first that 9 f m. In fTfl (2.3)] 
the following bound is proved: 

Aoo(g) > llog\x{Q)'' + 5Am^x{Q)\ - log A*. (16) 

o LZ 

The bound in (1161) holds for any rational point. 

If Q G -Eo(Qp) for all primes p then sum over all p, using (|T6l) 
and (|T5|) . to obtain 



kQ) > llog\x{QY + 54mMQ)\ = l\og\x{Q)\ + l\og\x{Qf + 54m% 

o o o 

using the product formula to write X]p<oo ^'-'S I^*Ip = 0- Since 

x{Qf > 27mV4, 

it follows that 

1 /27\ 111 f27 \ 1 
h{Q) > — log ( — j+— logm2+-logm^+-log ( — + 54j > -logm. 

If m ^ ±2 mod 9 then, according to [HI Page 180], 3P G ^o(Qp) for 
all primes p so we may put Q = 3P. The lower bound 

h{P) > ^^ogm 

results, which is stronger than the lower bound f[T^ . If m = ±2 mod 9 
but p has no prime factors congruent to 1 mod 6 then, by |I7|, Page 
180] again, we may take Q = 2P. This time we obtain the lower bound 

h{P) > — logm 

which, again, is stronger than the lower bound in f|T3|) . 
When 9 I m, we use [171 (2-4)], which gives 

Aoo(g) > I log \x{Q f + 6M'x{Q) I - ^ log A*, 

in place of (ITB]) . Assuming again that Q G -E'o(Qp) for all primes p and 
summing over p, we obtain 

HQ) > I fog KQY + 6M^x{Q) I = 1 log |x(Q) I + ^ log \x{Qf + 6M% 
8 8 8 

When 9 I m, x(g)3 > 3MV4 so 

HQ) > ^ fog (l) + ^ log + ^ log M' + l log + 6 
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which simphfies, upon setting M = m/9 to 

HQ) > ^ log m - ^ log 3 - ^ log 2. 

Again by [171 Page 180], we have 3P G i?o(Qp) for all primes p and 
putting Q = 3P yields 

KP) > -^^ogm - ^log3 - ;^log2 > 7^ logm - .0562 
27 3o 27 27 

Thus ( fT3l) holds when 9 | m. 

Finally consider the case when m = ±2 mod 9 but p has prime 

factors congruent to 1 mod 6. In this case Q = 3P G -E'o(Qp) for all 

primes p ^ 3. When p = 3, use [T7', (2.6)] to obtain, 



A3(g)>-^log3-^log|A*|3. 
Now sum over all p, using (ITSl) . (|T6l) and (|T7j) . to obtain 

/i(Q) > I log |a;(Q)^ + 54m2x(Q)| - ^ log 3. 
8 2 



(17) 



Since x(g)3 > 27mV4, 

1 /27\ 111 
/i(<5) > TTjlog ^ + — logm2 + -logm2 + -log 



24 



which simplifies to 



24 



27 

— + 54 
4 



2 logs. 



113 
/i(Q) > - log m - - log 2 - - log 3. 

Writing Q = 3P and dividing through by 9 gives the bound in (IT^ . □ 

4.2. Proof of Theorem 12.31 Assume that n is an index such that 
An has no primitive divisor; the following proof shows that n < 12. 
The proof comes in two steps. For the first step we will show that the 
assumption implies a divisibility statement of the following kind: 



An\2^'3^p{n)l[Ar_ 



(19) 



q\n 



where the product is taken over primes, and where p{n) denotes the 
product of the primes g > 3 which divide n, 

p(^) = n ^- 

3<g|7i 

This step is crucial in the sequel. A little parsimony here greatly re- 
duces the number of values m which need to be checked manually. 

To prove (fT9|) . let p be any prime dividing An- Assuming An has 
no primitive divisor, there is a term Ak with k < n such that p \ Ak- 
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Initially, suppose that p \ Q. If a is the rank of apparition of p in the 
sequence A then we know, by the discussion in Section [31 that p \ 
if and only if k = da for some d prime to 3. Since a < n, there is 
some prime q, necessarily distinct from 3, such that qa \ n, and hence 
p I An/q. Applying Lemma [3 ■2[ we have 

ordp(y4„) = ordp(y4|) + ordp(g) < ordp(y4a) + 1. 

Now consider the possibilities when p | 6, subdividing according to 
the cases in Lemma 14. 1\ beginning with the simplest. 
Case II 2 f m, 9 f m. 

For the case when p = 2 notice that there are no rational points on 
the minimal model with a; = mod 2. This is because the expression 
{27w? + l)/4 is odd and the equation y"^ + y = 1 has no solutions in F2. 
Hence the maximum value of ord2(v4„) is 2, using the transformation 
X = 4x, and indeed this is the 2-adic order of all the A^. In other 
words we may may take fi = in flTIJl) . unless n is prime, in which 
case /i = 2. 

Similarly, on the minimal model there are no points with x = mod 3. 
To see this, notice that (27m^ + l)/4 = —20 mod 81, if 3 | m, and the 
equation y"^ + y + 20 = mod 81 has no solutions. If 3 f m then 
(27m^ + l)/4 = 7 mod 81 but y"^ + y + 7 = mod 81 has no solutions. 
Again, using the transformation X = 4x, we see that 3 | An cannot 
hold and we may take A = in (fT9|) . 

The possibilities for the remaining cases can be summarized as fol- 
lows: 

Case I 2 | m, 9 \ m. If 4 f m then x = mod 2 does not hold while 
if 4 I m then a; = mod 4 does not hold. Using the transformation 
X = 4x, we see that 2 < ord2(A„) < 3 so we take /i = 1, unless n 
is prime, in which case /i = 3. For the prime 3, when 3 f m we find 
a; = mod 3 does not hold while if 3 | m we find x = mod 9 does not 
hold. Hence we may take A = unless n is prime, in which case A = 1. 

Case III 2 | m, 9 | m. If 4 f m then x = mod 2 does not hold 
while if 4 I m then x = mod 4 does not hold. Using the transformation 
X = 36x, we see that 2 < ord2(A„) < 3 so we take fi = I unless n is 
prime, in which case fi = 3. For the prime 3, we find a; = mod 3 does 
not hold. Hence we may take A = unless n is prime, in which case 
A = 2 (from X = 36x). 

Case IV 2 f m, 9 | m. We find x = mod 2 does not hold and 
neither does x = mod 3. Using the transformation X = 36x gives 
/i = A = unless n is prime, in which case /i = A = 2. 
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Now that ( |T9l) is established, take logarithms: 

log An < log p{n) + ^ log Ah + yu log 2 + A log 3, 

q\n 

where the sum runs over primes q. In each case, going to the minimal 
model means we must substitute An = v?an- This yields 

loga„ < logp(n) + ^logaH + 21ogu(tu(n)-l)+yulog2 + Alog3, (20) 

q\n 

where uj{n) denotes the number of distinct prime divisors of n. 
In order to apply Lemma (14.21) we will need the following bound: 

54M3 



log 



mxt 



< log9 = 21og3. 



(21) 



To prove (|2T1) note firstly that when 9 f m, we have m = M. Using 
either form of the minimal equation in Lemma 14.11 it follows that 
Xn ^ 27m^/4 and the bound in fl2Tl) follows at once. When 9 | m, 
M = m/9 and Lemma [4.11 gives > 3M^/4. Now (|2T1) comes out in 
exactly the same way. 

Multiplying ffTTl) by 2 gives 

2hn — - log 3 : log 



6 

and inserting (l2Ti) gives 



4 



1 + 



< log ar, 



2hn^ — - log 3 < loga„. 
Inserting this into (!20|) . we get 

2hn^ — -log3 < logp(r;,)+^ loga2i+2(co'(r;,) — 1) logu+/ilog2+A log3, 

q\n 

(22) 

On the other hand, multiplying (fT2|) by 2, replacing n by n/q, and 
inserting into ( l22l) gives 



2/m^ log 3 < logp(n) 

3 



Y^Uhf-) +^logM + 31og3j 

gin \ ^ ^ ' / 



+2{uj{n) -l)logM + plog2 + Alog3. 
Re-arranging gives 



2hn' 



1-^4 < logp(n) + ^uin) logM 
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+2{uj{n) - 1) logM + /ilog2 + ( A + 3uj{n) + - ] log 3. 



3 



To ease the notation, write 



q\n 
q prime 

Since M < m in each case, dividing by logm yields 



2/1^2 4 logp(n) +cu(n)log(27u2) +log(2'^3^+l/M2) 

-f{n) < -uj{n) + 



log m 3 log m 

(23) 

Case II 2 I m, 9 f m. 

In this case, u = 2 and fi = X = unless n is prime, in which 
case fi = 2 and A = 0. When n is composite, assume m > 40. The 
bound f|T3l) inserted into fl23l) now forces n < 12. Assuming m > 290, 
the bound f[T^ inserted into fl2^ forces n < 12. Thus we need to 
check all cases manually when m < 40 and all cases m < 290 when 
m = ±2 mod 9 and m has a prime divisor = 1 mod 6. This will follow 
in the next section. 

When n is prime, uj{n) = 1 but we need to take /i = 3 and A = 1. 
When m > 40, ([13]) and ([23D force n <7. When m > 290, ([HD and 
(1231) force n<7. 



Remark 4.5. To obtain these explicit upper bounds rigorously re- 
quires a little more than checking (125]) . because some of the functions 
that appear in (125]) are not monotonic. In fact (125]) implies a weaker 
inequality, where u{n) is replaced by its upper bound of log n/ log 2, 
p(n) is replaced by its upper bound of n, and /(n) is replaced by its 
lower bound of .547. The resulting inequahty yields a rigorous upper 
bound and the remaining cases can be checked manually. In none of 
the cases needed does this higher bound exceed 24 so the amount of 
extra checking is negligible. 

Now we summarize the arguments in the other three cases: 

Case I 2 | m, 9 f m. Here u = 2 and /i = 1 and A = unless n 

is prime, in which case yU = 3, A = 1. Applying Lemma 14.31 together 

with (123]) as before gives the same bounds as the previous case with 

the same manual checking to be done. 

Case III 2 | m, 9 | m. Here m = 6, /i = 1 and A = unless n is 

prime, in which 3 and A = 2. Inserting Lemma S]3] into (123]) . 

the same bounds arise. 
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Case IV 2 \ m, 9 | m. Here u = 6 and /i = A = unless n is prime, 
in which case /x = 2, A = 2. Inserting into (1251) and using Lemma [4.31 
gives exactly the same bounds as before. 

The proof of Theorem 12.31 is complete subject to checking various 
values of m; see Appendix lA.ll for the details. □ 



4.3. Proof that Z{W) < 14. Note that if a prime p is a primitive 
divisor of a term An, where n > 2, then it cannot divide with 
r < n. We have seen that a prime of bad reduction which divides An, 
with n > 1, also divides Ai, and so such a prime cannot be primitive 
divisor of An- If p > 5 is a primitive divisor of An then p is a prime 
of good reduction for E and so, by the results in Section [3l the rank 
of apparition of p in the sequence B is 3n. When p \ Q, the same 
conclusion follows by manual checking, using the same case-by-case 
analysis as in Section I4l2l 

The proof that Z(W) < 14 runs along almost the same lines as the 
proof that Z{A) < 12. Suppose p > 5 is a primitive divisor of An, and 
suppose that Then either p divides Wn, in which case it is a 

primitive divisor, or p"" \ 36mi?^ + C„, according to (jS]). By the identity 

(36m5^ + Cnf + (36m5^ - Cnf = m(6A„5„)^ 

which is simply ([T]), we have | SQmB^ — Cn, whence p"' \ 36mi?^. As 
An is prime to Bn, and as p > 5, we have p"" \ m. Thus Wn fails to 
have a primitive divisor just in case the 'primitive part' of An divides 
m, or 

An\m2''3^p{n)Y[A!i. (24) 

q\n 

The arguments of the previous section now apply mutatis mutandis 
to show that 

2hn'^ 4 logp(n) +cj(n)log(27M2) + log(2'^3^+§/M2) 

/(^) ^ 1 + TT^l'^) "I \ • 

log m 3 log m 

(25) 

When n is composite, assume m > 40. The bound f[T^ inserted 
into f l25l) now forces n < 14 in all cases. Assuming m > 290, the 
bound (IT^ inserted into (125!) forces n < 14. This completes the proof 
that Z{W) < 14, subject to checking various values of m. The values of 
m required to complete the argument are checked in Appendix lA.li □ 
Thus the proof of Theorem 11.11 follows once we establish that Wn 
has a primitive divisor for all 2 < n < 14. To that end we now turn. 
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5. Primitive divisors in specific terms 

In the present section we concern ourselves with the question of, for 
fixed n, which sequences W fail to have a primitive divisor in the nth 
term. For many small values of n we can show that there are no such 
sequences, and this will bridge the gap between the work in Section H] 
and the goal of Theorem ll.il To ease notation, we will frequently write 
D in place of — 432m^. 

Proposition 5.1. Let W be a sequence as defined above. Then Wn 
has a primitive divisor for each 2 < n < 14. 

The proof will come in several pieces. For n > 5 not divisible by 3, 
the methods developed in [14] suffice to treat this problem. Although 
we make use of the special form of = — 432m^ , much of the argument 
will work for general Mordell curves. Indeed, we reduce the proposition 
to the checking of finitely many elliptic divisibility sequences arising 
from Mordell curves, none of which turn out to be curves of the special 
form under consideration. When n > 6 is divisible by 3 some problems 
arise, but only a slight modification of the method of [H] is needed. 
The cases n = 4, n = 3, and n = 2, treated in Section El are dispatched 
largely through ad hoc means, although the spirit of the proof remains 
the same. 

5.1. Division polynomials. We exploit, as above, the short Weier- 
strass form of the equation, and will in fact show that Wn has a prim- 
itive divisor coprime to 6m for the n listed in the lemma. We will 
consider the curves C and E, as defined in ([1]) and (jl]) respectively, 
bi-rationally equivalent by the map defined in iQ. There are (see, for 
example, |i24j) rational maps (pn.ipm^n £ Q.i.E), the function field of 
E., such that, for each G Z and Q G -E'(Q), 

In Section |H we exploited the existence of primitive divisors in the 
sequence A. Here we use the properties of the division polynomial 
tpn, much as in |l3j, to establish the existence of primitive divisors in 
Wn- An approach employing the polynomials 0„ would be similarly 
successful, and would be in keeping with the fiavour of Section HI but 
such an approach also encounters serious computational difficulties. 

We will show that for all n under consideration (save those treated 
later as special cases), there is a prime p dividing _B„ such that p \ GAkBk 
for any k < n. Such a prime will divide Wn, as gcd(i?„, C„) = 1, but 
not Wk for k < n. That is, p will be a primitive divisor of Wn- It turns 
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out that this amounts to showing that there is some prime p \ ^{E) 
such that ordp(F„(Q)) > 0, for a certain F„ G Q(-E) constructed below. 
The hypothesis that there is no such prime will lead us to a non-trivial 
solution of a certain Thue-Mahler equation which depends on n but, 
critically, not on Q or m. By explicitly solving the Thue-Mahler equa- 
tions in question, we will treat all bar three cases of Proposition I5.1[ 

As Q{E) = Q{x,y), with = + D (for fixed D = -432m^), 
we may write 0„, ipn, and Un as polynomials in x and y, where y 
occurs only to the first power. In fact (see [24j), we may write 0„ as 
a polynomial in x alone, and either ipn or y~^ipn similarly, as n is odd 
or even respectively. Note that, in all cases, ip"^ may be written as a 
polynomial in x, and we will view (pn and ip'^ as elements of Q[x]. In 
fact, if E[n] denotes the kernel of multiplication by n in E{Q), 



Note that each linear term on the right side occurs precisely twice, 
except those corresponding to T G -E[2]. 

As our proof relies on the properties of ipn, it is incumbent upon us 
to make a few observations (these remarks are made as well in [Hj). 
Note that if 3 f n, then {x, y) G E[n] implies both x 7^ and 



where C is a primitive 3rd root of unity. In particular, ipn or 'ipn/y, 
depending on the parity of n, is a polynomial in x^; ip"^ G Q[x^] C Q[x]. 
When 3 I n, on the other hand, ip"^ G x^Q[x'^]. 

Although ipn depends on m, our stated aim is to construct, from 
each term Wn without a primitive divisor, a solution to a Thue-Mahler 
equation that is independent of m. Note that the points ony'^ = x^ + D 
and those on y"^ = x^ + D' are related by the scaling map 



Thus, the dependence of fl26l) on D is transparent. Along the same 
lines as jH], we observe that ip'^ may be written, over Z, as a binary 
form in x^ and AD. We will abuse notation somewhat, and denote this 
form by ip"^ as well, so that 



In general, ip'^ G Q(-E') is not irreducible. Aside from being a square 
when n is odd, the function has several obvious factors. Let CooiT) 



(26) 



T(^E[n] 



{x,y), {Cx,y), {C'^x,y) G E[n\ 




i^l{Q)=^Pl{x'AD). 
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denote the order of the torsion point T G E{Q), and let 

F^{Q)=e'{n) n {x{Q)~x{T)), (27) 

TeE{Q) 
e^{T)=n 

where 

\p if 77, = is a prime power 
e[n) = < 

I 1 otherwise. 

It is clear from (l26l) that 

d\n 

Note that, when n ^ 2, each term on the right of ( 1271) occurs precisely 
twice, allowing us to define Fn G Q[x] implicitly in this way. 

The functions F„ may be viewed as the elliptic analogues of the 
cyclotomic polynomials exploited in [3l O [21] to treat the analogous 
problem for Lucas sequences. By the same arguments as in the case 
of ipn, the product in (1271) defines a binary form over Z in and 4D, 
at least when n ^ 3. We will, then, write F„(x^,4D) for the product 
above. We will frequently pass between the rational function Fn{Q) and 
the binary form Fn{Al,4:DBf). It is worth noting that, as Ai and Bi 
are co-prime, the primes appearing to a positive power in Fn{Q) occur 
to the same power in F„(y4f , 4Di?^), except possibly those dividing 
e{n). 

5.2. Division polynomials in finite fields. Note that if p is a prime 
of good reduction for E, then the equations from [21] define the divi- 
sion polynomials in the same way for E{¥p). In particular, if ep{Q) is 
the order of the image of the point Q G -E'(Q) in the group -E'(Fp), then 
ep{Q) is precisely the rank of apparition of the prime p in the sequence 
(V'n(<5))n>i, that is, the smallest n such that ordp(^/'„(Q)) > 0. Equiv- 
alently, ep{Q) is the unique n such that ordp(F„(Q)) > 0. Note, on the 
other hand, that x{Q) = mod p implies that Q is a point of order 
3 in E{¥p). More generally, then, if n is the rank of apparition of p 
in the sequence (0n(Q))n>i, it follows that ep{Q) = 3n. The converse 
may fail, of course, as there may be points of order 3 in E{¥p) other 
than those with x = 0. The following observation will be useful. 

Lemma 5.2. Let Q G E{Q), let An and En he defined as above, and 
let p \ A{E) . Then the following hold: 

(a) if ordpiipniQ)) > 0, then p \ B^; 

(b) if ordp{(f)n{Q)) > 0, then p \ A^; 
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(c) if OTdp{Fn{Q)) > and 3 \ n, then p \ Bn and p \ GAkBk for 
k < n. 

Proof. Note that results similar to (a) and (b) are derived in [T|[T^[T6]. 
In particular, the result in [1] is stronger: the hypothesis is merely that 
the reduction of Q modulo p is not singular. The lemma is demon- 
strated here for completeness. 
As 

An ^ MQ) 

the conclusions of (a) and (b) can only fail if ordp((/)„(Q)) and ordp (■?/'„ (Q)) 
are simultaneously positive. If this is the case, then ordp(?/'fc(Q)) > 
for = n lb 1 as, by definition (see pi]). 

By an easy modification of Lemma 4.1 of |29j (for details on this mod- 
ification see [16j), we have ordp(?/'3(Q)) > and oidp{ipi{Q)) > 0. Com- 
puting the resultants of the binary forms ip3{x^, 4D) and ^/'|(x^, 4Z}) in 
Z[D], we see that this can happen only if p \ 6m. That is, if p \ A{E). 

Now suppose OTdp{Fn{Q)) > 0. Part (a) ensures that p \ Bn, and 
the image of Q in E{¥p) has order exactly n. If p | B^ for k < n, then 
OYdp{ipk{Q)) > 0, and hence ep{Q) < k < n, a. contradiction. Similarly, 
if p \ Ak then we have ordp(0fc(Q)) > 0. Hence the prime p has finite 
rank of apparition in the sequence {4>n{Q))k>o, and so n = ep{Q) is 
divisible by 3. We assumed it was not. □ 

5.3. Thue-Mahler equations. The factors F„ are the binary forms 
exploited in [Til [TB] to show the existence of primitive divisors in spe- 
cific terms of elliptic divisibility sequences arising from elliptic curves 
in short Weierstrass form. We will use a similar approach here to 
construct primitive divisors of the sequence B. The following, simple 
observation will be used repeatedly. 

Claim 5.3. Let 

s = Al/ gcd{Al AD), t = ADBl/ gcd{A\, AD). 
Then for all primes p > 5 dividing gcd(Ai, D), we have 

ordp(s) > ordp(t) = 0. 
Furthermore, for all primes p > 5, the quantity ordp(t) is even. 
Proof of the Claim. Note that, as m is cube-free, we have 



ordp(m) G {1,2} 
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for all p > 5 dividing m. Suppose that p \ Ai (and hence p \ Bi). If 
ordp(m) = 1, then 

OTdp{A\) > 3 > 2 = oidp{ADBl) = ordp(-432m2). 

Thus ordp(s) > ordp(t). 

If, on the other hand, ordp(m) = 2, then divides the right-hand- 
side of 

= Al- A32m'^Bl 

But then divides Af — A32m'^Bf, as p divides the left-hand-side of 
the above to an even power. Now p'^ \ Ai, hence 

ordp{Al) > 6 > 4 = ordp(-432m2), 

and again ordp(s) > ordp(t). The fact that ordp(t) = simply follows 
from s and t being, by construction, relatively prime. 

Now suppose that p > 5, and further that ordp(t) 7^ 0. We have just 
shown that this ensures p\ Ai, and so 

ordp(t) = ordp(— 432m^) = 2ordp(m) = mod 2. 

□ 

We are now in a position to present the main tool in the proof of 
Proposition 15.11 for the values of n > 5 which are not divisible by 3. 
The values n = 2, 4, and those values divisible by 3, require a slightly 
more careful treatment. 

Lemma 5.4. Suppose that Wn has no primitive divisor, for n > 5 not 
divisible by 3. Then, for s and t defined as above, Fn{s,t) is divisible 
only by primes dividing 6e{n). 

Proof. We have seen that if p f ^{E) and ordp(F„((5)) > 0, then p is a 
primitive divisor of Wn- Thus, if Wn has no primitive divisor, we must 
have OYdp{Fn{Q)) > only for p | A{E). Clearing denominators, we 
see that the only primes dividing the integer Fn{Al,4:DBf) are those 
of bad reduction for E or, possibly, those dividing e(n). If p f 6e(n) is 
one such prime, then p \ D. But then 

Fn{Al,ADBf) = e(n)A?'^'^^" mod p. 

In particular, we can only have p | Fn{Al,4:DBf) if p | Ai. In the 
latter case we have p | gcd{Ai, D) whence, by Claim [375| we obtain 
ordp(s) > ordp(t) = 0. Note that -F„(0, 1) = ±1 as in [H] (see also the 
tables in Section 1X121) . Thus 



Fnis,t) = it^egK ^0 mod p. 
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In particular, the prime divisors of Fn{s,t) are at most those of 6e(n): 

□ 

Note, for the purposes of impending computations, that we can say 
somewhat more. If n = p° is a prime power (and so e{n) = p), then 
OYdp{Fn{Q)) > and p \ A{E) together imply p \ m. In this case 
D = mod p^, ergo 

FniAlADBf) = pAl'''^'-^"^ mod p\ 

Thus if p f Ai, we have ordp(-F„(A'^, ADBf)) = 1. If, on the other hand, 
p I Ai, we have ordp(s) > ordp(t), and so j9 f t). In the notation 
of the proof, then, we may take 7 G {0, 1}. We will see, by examining 
the individual forms, that we can control the exponents of 2 and 3 by 
elementary means as well. 

5.4. Solving the Thue-Mahler equations. We have reduced the 
proof of Proposition 15. II to treating the special cases n = 2, 4 and 3 | n, 
as well as solving a number of Thue-Mahler equations. Although the 
proof appears to require us to find all solutions to the Thue-Mahler 
equations 

F„(s,t)=2"3^e(nr, 

we have already seen that we need only consider equations wherein 
7 e {0, 1}. Although we can, a priori, restrict the exponents a and (3 
by employing such techniques as lower bounds on linear forms in p-adic 
logarithms, it turns out that we may also do so by more elementary 
means. 

For the various n under consideration, consider F„ mod 2, and note 
that 

Fn{s,t) = 1 mod 2 

whenever gcd(s, t) = 1 (note that these forms are available below for 
examination). Similarly, we may consider F„ mod 3 and reduce the 
possible values of P to three: 0, deg(F„), or |deg(F„). Thus, our 
Thue-Mahler equations may be reduced to 6 or 12 Thue equations 
(depending on whether or not n is a prime power) of the form 

F„(s,t) = (-l)^3^e(n)^ 

with 5,7 e {0,1} and f3 G {0, deg | deg(F„)}. 

We may also disregard several possible solutions in advance. For ex- 
ample, some of these equations possess a solution (s, t) with t = 0. Such 
a solution cannot arise from a pair (s, t) as constructed in Claim 15. 3[ 
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however, as this would necessitate either m = or i?i = 0. Further- 
more, we have from Claim that ordp(t) is even for all p > 5, and so 
we may ignore solutions (s, t) that fail to have this property. As 

tgcd{AlAD) = -A32m^Bf < 

and 

(4s + t) gcd(A?, AD) = A{Al + DBf) = ACf > 0, 

we may conclude as well that t < and 4s + 1 > (here, notice that 
Ci = only for points Q of order 2, which are not the type under 
consideration). Any solution not satisfying these inequalities may be 
discarded as well. Solutions wherein s = correspond to points Q of 
order 3, and the solution (s, t) = (1, —1) gives rise to the point (12, 36), 
a point of order 3 on = — 432. We shall call solutions falling into 
the above categories expected. 

A computation in PARI/GP [IS] shows that there are no solutions to 
any of the Thue equations above, other than (possibly) these expected 
solutions. In the appendix, we list the binary forms Fn for the various 
values of n, in order that the reader may confirm these findings. We 
should note that the Thue equation solver in PARI/GP assumes, by 
default, the truth of the Generalised Riemann Hypothesis. This default 
was overridden, and our results verified unconditionally. The most 
strenuous computation arises in the case n = 11, in which a Thue- 
Mahler equation of degree 20 must be treated. (The binary form arising 
in the case n = 13 is of greater degree but factors, and so the equations 
above may be treated by elementary means.) Even this computation, 
however, took well less than a minute (on a 1.83 GHz MacBook with 
512MB of RAM). 

5.5. The case 3 | n, but n ^ 3. Some care must be taken when 3 | n 
and n ^ 3, but the methods are not fundamentally different. 

If Fn G Q{E) is defined as above, where n = 3k, and if p is a 
primitive divisor of Ak, then we have ordp(F„(Q)) > 0. It is not clear, 
then, that the primitive divisor of Bn which we prove to exist by the 
method above fails to divide Aj.. We must modify our argument to 
show that Bn has a primitive divisor distinct from those coming from 
Ak. The binary form F„ factors in this case, and we will show that it 
suffices to show that one of the factors has a prime divisor other than 
2 or 3. 

Let H (1E{Q) be the group 

|(0, /D), (0, -Vd), o} ^ Z/3Z. 
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Note that for D = -432m'^, H C E{Q{^/^)), although much of what 
is written here apphes for all D. For Q G E{Q), let e{Q, H) denote the 
least k such that kQ G H, if one exists, and e{Q, H) = oo otherwise. 
Set, for n divisible by 3, 

Gi{Q)= n i<Q)-<T)). 

TeE{Q) 
e{T,H)=n/3 

As H is fixed by the action of Galois on Q, we have immediately that 
Gn G Q(-E). If Gn vanishes at Q, then An = 0, and so we see that 
Gn I (pn- If we define Hp C E(¥p) to be the analogous subgroup (for 
p \ A{E)), and ep{Q,H) to be the analogous value, then ep{Q,H) 
is the rank of apparition of p in the sequence A, as discussed above. 
Arguments similar to those before show that Gn may be written as a 
binary form in and 4D with coefficients in Z, and we will denote 
this form by Gn{X, Y). For example, 

Gg{X, Y)=X^ - 2AX^Y + SXF^ + Y^. 

Note that the roots of Gn are points on E{Q) of order precisely n. Thus 
Gn divides Fn. Let 

Fn{X,Y)=Gn{X,Y)Fn{X,Y). 

Lemma 5.5. Suppose that Wn has no primitive divisor, for 3 | n. 
Then Fn{s,t) is a {2,3}-unit. 

Proof. Suppose that p \ A{E) and ordp(F„((5)) > 0. It follows that 
OTdp{Fn{Q)) > 0, and so, just as in Lemma [521 we have p \ Bn and 
consequently p \ Wn- We have, by the same argument as is Lemma [5^ 
that p \ Bk ioT k < n. Suppose p \ Ak, for some k < n, and suppose 
without loss of generality that k is the least such value. Then we have 
^p{Q) = 3fc, and so n = 3k. We have as well that ep{Q,H) = n/3, so 
OTdp{Gn{Q)) > 0. But Fn and Gn can have no common roots modulo 
p. If Q were such a root, then x{nQ) = mod p while also x{nQ)^ = 
—4D mod p, clearly impossible if p\ D. So we must not have p \ A^. 

If, on the other hand, p \ A{E) and p > 5, we have p | m. As above, 
we have ordp(s) > ordp(t) = and so, as Fn{l, x) is monic, p f F„(s, t). 
The result is proved. □ 

Note that the lemma is true in the case n = 3. The binary form F^ 
has degree 1, however, and so we cannot proceed in the same way as we 
will for n > 6. As in the case of Lemma [5.41 some care must be taken 
with the exponents. We note that Fn{s,t) = 1 mod 2 for all relatively 
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prime s and t, while 

{0, 3, or 5 if n = 6 

0,9, or 13 if n = 9 

0, 12, or 18 ifn = 12. 

Computations in PARI/GP are as in the previous case, and reveal 
no unexpected solutions. 

6. Three special cases 

We have, thus far, shown that Wn has a primitive divisor for all 
n > 5. We treat the cases 2 < n < 4 here, completing the proof of 
Theorem II. 1[ 

6.1. The cases n = 3 and n = 4. To treat the case n = 3, we 
show that ^3 must have a prime divisor not dividing 6mAiA2. Such 
a divisor, being a prime of good reduction for E, must divide W3, but 
cannot divide Wi or W2. Suppose, to the contrary, that A3 has no 
prime divisor other than those dividing 6mAiA2. Comparing (^3 to '02 
and 02, we see that then any prime p with ordp(03((5)) > must be a 
prime of bad reduction for E. If p > 5 is a divisor of m we see, just as 
above, that 

03(s, t) = - 2AsH + ?>st^ + 

is not divisible by p. So ^3(5, t) is a {2, 3}-unit. Exactly as in previous 
cases, we reduce this to 

- 2AsH + ?,sf + = ±3^^, 

with P G {0,3,4}, by considering the possible values of 03 modulo 2 
and 3. Computation produces no unexpected solutions. 

For the case n = 4, we will show that ^4 has a prime divisor not 
dividing 6mAiA2A3. As this divisor is a prime of good reduction for E, 
it divides W4, witnessing that W4 has a primitive divisor. We consider 
one particular factor of (j)l{X,Y). Let 

F;(X, Y)=X^- 134X^y - SAX'^Y^ - 32XY^ - 2Y^. 

One can verify that this binary form divides 04, and (by the com- 
putation of some resultants) that F^{Al,4:DBf) has no prime factors 
in common with Ai, A2, A^, Bi, B2, B^, save possibly some divisors of 
A{E). By a final application of Claim [531 we see that Fl{s,t) is a 
{2, 3}-unit. Solving the implied Thue-Mahler equations, one finds no 
unexpected solutions. 
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6.2. The case n = 2. It remains to check that, for any m and Q, the 
term W2 has a primitive divisor. From (P) we have, in the notation of 
previous sections, that 

U2_ ^ 36m^l{Q) + uj2{Q) 
W2 ~ 6MQ)MQ) ' 

where 

(f)2 = x{x^ -8D), '^l = 4:{x^ + D). 

We suppose that every prime dividing W2 also divides Wi, and hence 
6A1B1. If p > 5 is any prime with 

oidpiMQ)) = ordp(4A? + ADB^,) > 0, 

then p is a primitive divisor of W2, or p \ Qm. Assuming that the former 
is not the case, then, we see that (4s + 1) is a {2, 3}-unit, employing 
Claim [573] as in the cases above. That is, we know that p \ (AAf+ADBf) 
only if p I gcd(Ai,m), in which case ordp(s) > ordp(t). In fact, we can 
say slightly more. Just as in Section 15. 4[ we may conclude from the 
inequality 

{AAl + ADBl) = ACf > 

that (4s + t) is a positive {2, 3}-unit. 

Now consider the primes dividing v4f- 81)5^. Just as in the previous 
case, if p > 5 is a prime divisor of this expression, and p is not a 
primitive divisor of W2, then p \ m. If this is the case, then another 
application of Claim [5731 tells us that (s — 2t) is a {2,3}-unit. Solving 
two linear equations yields 

gs = Wi + 2W2 

gt = -AWi + W2, (28) 

where Wi and e Z are {2, 3}-units, and W2 > 0. It is in which 
we are most interested. The second part of Claim 15.31 tells us that 
ordp(t) is even for any p > 5. Specifically, then, t = —dx'^, for some 
positive d \ 6 and some x G Z (recall that t < as per Section r5.4p . 
Equation fl28|) is now the representation of — dx^ as a sum or difference 
of {2,3}-units. We also know that gcd{Wi,W2) \ 9, as gcd(s,t) = 1. 
We may use this information to solve the above system of equations 
for all possible values of s and t. First, a lemma. 

Lemma 6.1. The only integral equations a = b + c such that 

(a) c < 0; 

(b) b and c are {2,3}-units; 

(c) gcd(6, c) = 1; and 

(d) a, 2a, 3a, or 6a is a perfect square 
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are the following: 



1^ 


= 2^ 


-1, 


t 


= 3^ 


-2^ 


12 


= 2^ 




2.1^ 


= 3^ 


-1, 


12 


= 3^ 




2.2^ 


= 32 


-1, 


12 


= 32 


-2^ 


2.11^ 


= 3^ 


-1, 


52 


= 33 


-2\ 


3.1^ 


= 22 


- 1. 



Proof. Supposing that a = b + c is one such equation, we may multiply 
both sides by sufficient powers of 2 and 3 to obtain an equation of the 
form 

g2 = ^3 _ 2^3^^^ 

where r is a {2,3}-unit. Thus, if we write ^ = fiQ + 6/^1, where 
< /^o < 6, and write u = uq + Qvi similarly, we see that 

r q 



22^*1321/1 ' 23Ati 331/1 

is a {2, 3}-integral point on the elliptic curve = —2'^°3'^'\ with the 
additional property that X is a {2, 3}-unit. Using MAGMA |J8J to find 
all {2, 3}-integral points on each curve of this form (for < fio, ^0 < G), 
and tracing these points back to the original equations, we have our 
result. □ 

It is now a simple matter to solve the above system of equations for 
all possible s and t. If gcd(iyi, W2) = 9, then 

-t = - 

is an equation as in the lemma. Moreover, as ^Wi is divisible by 4 and 
W2 > 0, most of the equations listed in the lemma may be disqualified 
immediately. We are left with the possibilities (VTi, 1^2) = (9,27) or 
(9,9). These yield (s,t) = (7, —1) or (3, —3), the latter of which may 
be discarded as gcd(s, t) = 1 by construction. 
If gcd{Wi,W2) = 3, then 

-3t = -Wi - -W2 
3 3 

defines an equation as in the lemma, yielding again the solution (3, —3) 
for (s,t). The case gcd(Vri, W2) = 1 yields no new solution. 

The only solution, then, that might contradict our claim, is that 
corresponding to (s,t) = (7, —1). One may trace this solution back to 
the point (2, —1) on the curve + = 7. While the terms Af + DBf 
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and Al — 8DBf do both turn out to be {2, 3}-units in this example, we 
may simply compute Wi = 1 and W2 = 3 to see that W2 does, indeed, 
have a primitive divisor. 

This completes the proof of Proposition 15.11 and hence the proof of 
Theorem II. 1[ □ 

Appendix A. Computations 
In this section we discuss some of the particulars of the computations. 

A.l. Small values of m. As required in Sections 14.21 and 14.31 we 
computed manually the Zsigmondy bounds Z{A) and Z{W) for all the 
cases needed with m < 290. From the analysis given in Section 14.21 
it is sufficient to consider all m < 40 (in fact we go out to 50) but 
only those 40 < m < 290, for which m = ±2 mod 9 and which possess 
a prime divisor congruent to 1 mod 6. In all cases when the rank is 
positive the rational torsion group is trivial. 

For rank-1 curves, the values of Z{W) are given in Figure [H The 
values of Z{A) are much simpler: Z{A) = except in the case m = 7, 
when Z{A) = 2. 

In rank-2 the situation is slightly more complicated. However, pro- 
vided h{Q) > 0.1, equation fl2^ yields the Zsigmondy bound Z{A) < 12, 
and equation (!25|) yields the Zsigmondy bound Z{W) < 14. In all 
cases, there were no non-trivial rational points with h{Q) < 0.1 so 
no further checking was necessary. In Figure [2] we list generators for 
each of the rank-2 curves in the range. These were looked up in Cre- 
mona's tables [7], when the conductor was below 10^, or computed 
using MAGMA [TS] for larger conductors. 

A. 2. Division Polynomials. Various specific binary forms are used 
for computations in Section [5l and these are reproduced in Figures [3] 
and m To save space, only the coefficients have been recorded. The 
line 



F [vd,Vd-i, ■■■ ,Vi,Vo 



in the table is to be interpreted as the statement 

d 

F{X,Y) = J2'ViX'Y'^~'. 

Note that, as a result of the action of complex multiplication on E, the 
polynomial ipp factors when p = 1 mod 3. In actual computations, this 
fact can be exploited to great advantage, allowing the corresponding 
Thue equations to be solved by elementary means. 
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Figure 1. Zsigmondy bounds in rank-1 
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7 


[84,756] 


1 


9 


[36,108] 


1 


12 


[52,280] 





13 


[52,260] 





15 


[49,143] 





17 


[84,684] 





20 


[84,648] 





22 


[553/9,4085/27] 





26 


[156,1872] 


1 


28 


[84,504] 


1 


31 


[217,3131] 





33 


[97,665] 





34 


[2733,4455] 





35 


[84,252] 


1 


42 


[172,280] 





43 


[129,1161] 





49 


[196,2548] 





50 


[8148/27,138736/27] 





51 


[5473/36,333935/216] 





58 


[9444/27,173600/27] 





61 


[732,19764] 
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79 


[316,5372] 





97 


[388,7372] 





133 


[228,2052] 


1 


151 


[4228/9,261532/27] 





169 


[2028,91260] 


1 


223 


[1561/4,49283/8] 





241 


[6748,554300] 


1 


259 


[777,20979] 





277 


[5817/1^11261/8] 





286 


[58812960] 


1 



FlGL'llE 2. Rank-2 generators 
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[156,1908], [228,3420] 


37 


[84,36], [148,1628] 


65 


[129,567], [156,1404] 


124 


[372,6696], [2356,114328] 


182 


[273,245 7], [364,5824] 


209 


[1596,63612], [532,11476] 


218 


[1308,47088], [13881/25,1534221/125] 


254 


[16257/ L2072385/8], [508, 10160] 
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Fs I [5, 95, -15, -25, -1] 

F7 [7, 986, -2681, -12964, 3626, -1519, -686, -49, 1] 

Fs [2, 616, -7336, -1544, -3430, -4124, -952, -104, -1] 

Fio [1, 1173, -55284,29380,-368055,-1404072,-862941, 542232, 

. . . -104805,-7070,-474,-177,1] 
'11 [ 11, 23221, -1153603, -62045313, 66133914, -1596123771, 

. . .-8579472693, -4760052033, -22319781, 8054721004, 
...10595519759, 4869514969, 1106263389, 189881835, 
. . . 59389374, 17393277, 2270301, 102729,605,-242, -1] 
^13 [13, 74737, -10304874, -1459820466, 7383882519, -294761888811, 
. . .-3649379851026, -327751614216, 3634612800273, 
...75587434125411, 206422282971957, 165623202699903, 
. . . 77423927253309, 50317031121903, 70684315657137, 
. .. 64207462488471, 30461492791431, 8167061938581, 
. . . 1237534488021, 33446767107, -47530886481, -16133119236, 
. . .-2480541102,-183218139, -6445998, -217503, 

...-22815,-338, 1] 
14 [1, 8826, -3182349, 27544616, -1267563423,-29876807793, 

. . .-73452197357, -534368475927, -321414204609, 
. . .-159623734993, -250499094747, -930524257131, 
. . .-1172171589176, -509647490898, -20486729571, 
...61406271479, 22270327506, 3403598121, 263510632, 
. . . 15278739, 2663808, 488510, 19851, 537, 1] 

Figure 3. The binary forms F„ 



F12 



[1, 57, 3, 1] 

[1, 657, 6111, -3318, 19647, 12033, 3972, 684, 9, 1] 

[1, 3630, -28608, 392908, 212553, 1121508,168108, 62712, 

. . . 69507, 32782, 3684, 12, 1] 



Figure 4. The binary forms F„ 



